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Abstract 



Continuing our work , where a exphcit formula for the two-point functions of 
the two dimensional Z-invariant Ising model were found. I obtain here different 
results for the higher correlation functions and several consistency checks are 
done. 



1 Introduction 



In this paper I want to apply the the vertex operator method (VOM) for the calculation 
of different correlations functions of the two dimensional Z-invariant Ising model (IM). Z- 
invariance is the property that allows to introduce inhomogeneities in the lattice in conform 
with the integrability of the model. In particular when we put all inhomogeneities equal zero 
we obtain a homogeneous Ising model in a square lattice. Also by choosing appropriately 
the inhomogeneity parameters its is possible to obtain other regular lattice configurations, 
like a triangle, checkerboard and others (see for examples). The VOM was introduced 
by the Kyoto group in the case of the six-vertex model [0 and then extended to a variety 
of other integrable lattice models. The novelty of this approach is the treatment from the 
beginning of the lattice model in the thermodynamic limit, these is made possible because 
of two fundamental concepts: the corner transfer matrix (CTM) of Baxter and the lattice 
vertex operator (VO) introduced by the Kyoto group, see for a review. Whereas the 
CTM allows the calculation of the one-point functions the VOs permits in connection with 
the CTM the calculation of higher correlation functions (CF). Physically they correspond 
to inserting a semi- infinite inhomogeneity line at bulk of the lattice. This inhomogeneity 
is described by a spectral parameters that enters the parameterization of the Boltzmann's 
weights of the model such that they satisfy some integrability conditions (in the case of 
the IM they are called the star-triangle relations), which imply that the diagonal transfer 
matrices of the model commute and also allow the possibility of choosing different spectral 
parameters in the definition of the usual transfer matrices p|. 

In the paper the author obtained the two-point functions in the general Z-invariant IM 
whose vertex operator algebra was developed in . These were given in forms of pfaffians of 
well known elliptic functions over the inhomogeneity parameters. Because of its importance 
for the calculation of the higher correlation functions I review this result in Sect. |^ (see also 
Sect. [4.1D and in the Appendix I give explicit details of the calculations. After reviewing 
some important definitions in Sect. 0, then in Sect. 3^ I give a general integral formula 
for all higher correlations with an even number of spins. Particular cases of this formula 
are treated next. Section |^ contains our principal results, the first part review the general 
two-point functions as already mentioned, in Sect. ^]2| the homogeneous CF with adjacent 
diagonal spins are found in the form of generalized Toeplitz determinants and in Sect. [4.3| I 
find general homogeneous diagonal four-point functions, obtaining a general determinantal 
expression. In Sect. |^ I consider finally some particular cases of three-point functions. A 
conclusion section ends the paper. 



2 Fundamental definitions 
2.1 Elliptic parameterization 

Let us begin by introducing a number of definitions for the Z-invariant IM; the model is 
defined initially on a square lattice rotated by 45 degrees, the spin variables taking values 
or —1 live on the vertices of the lattice the coupling being between nearest neighbors. 
Assume first two different interaction energies (measured in kT units) with > in the 
NW-SO direction and L > in the NO-SW direction between two spins as in the figure |^ 
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Figure 1: Boltzmann Weights for the IM 



In the figure the two hues indicate one horizontal uh and one vertical uy spectral param- 
eters, their difference u = uh — uy ''is, used to code the couplings energies under the critical 
temperature by mean of Jacobian elliptic functions, 



where sum is a Jacobian elliptic function with half-periods /, iJ', k is the modulus, < 
u < I' (see e.g. and Glashier's notation for the Jacobian elliptic functions is used, 

nsM := l/snti, dsM := dnti/snw. These parameterization already introduced by Onsager is 
a consequence of the integrability conditions, in the case of the IM the so called star-triangle 
relations, which imply as well known, that diagonal matrices commute when the difference 
of the spectral parameters from site to site remain constant (see pag. 374 in Baxter's book). 
The module k = (sinh(2i^') sinh(2L))~^ does not depend on the spectral parameter and 
tends to 1 when the systems goes critical. It is possible to write a similar parameterization 
for temperatures higher than the critical but it will not be necessary due to the duality 
between order and disorder operators to be introduced below, which imply that the disorder 
CF under the critical temperature correspond to the order CF above it. 

2.2 Vertex Operator Algebra 

I review now the fundamental results of the paper where the vertex operator algebra for 
the IM was found. Dividing the diagonal lattice in four equal quadrants to apply Baxter's 
CTM method we have two cases according to wether the central point belongs or not to the 
lattice; this defines two CTM's Fig.H that act on the two sectors of the model, NS and R, 
respectively (in analogy to conformal field theory). Diagonalising the adjoint action of the 
CTM Hamiltonians we can describe the two sectors as Fock spaces (called Ti^^ and Ti^) of 
two irreducible representations of free fermions (apart from a different normalization). We 
have: 



sinh(2is:) = -i; 
sinh(2L) = ik~ 



isniu, cosh.{2K) = cniu, 
~^nsiu, cosh(2L) = ik~^dsiu, 



(1) 




are the respective 
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There are two vertex operators (VOs in the following) depending on which sector they 
intertwine: 




Figure 2: The CTMs for the IM 

The Vertex Operators (VO) satisfy the following elegant intertwining relations with the 
fermions 



a^.K) = ^f(C)$.K). (3) 
Here we use the simplified notation of |[ $^^'"(0 = $'"(0, 

a being equal to +1 or — 1 

distinguishes the two sub-sectors (with even or odd parity) of the NS sector and $NScr(C) = 
$o-(C), also we define [2;] := f{z) = —^/ksnv with z = exp{i7Tv/I) and we have introduced 
the generating functions 

= E A^) = E<^^"^ (4) 

rez+i sez 

and 
where 

[^]ns = y^cnt; = E 

rez+i 
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■kJ' 

rjr = X +x , x = exp(— ^yj. 

By crossing symmetry i.e. u u — I' oi in the multiplicative language — ^ C "we can 
obtain (using also the transformation f{x'^z) = l/f{z)) similar relations with the parameter 
of the VO replaced with ( but in the calculation of the CF in the next sections I will need 
only the expressions above. Finally the VOs satisfy the following unitary relations: 

5^$.K)$'^(C) = ^7^xid«H,, (6) 

<l>'^K)<f.(C)=^?'''xid„Ns,., (7) 

where the scalars and g^^ can be given as infinite products in the variable x but as we 
will not need them for the calculation of the CF we leave them undefined, similarly the 
expectation values of the VOs in the respective fermion basis are known but will not be 
needed (see for their explicit values). 



3 Even correlation functions 



In this section I will give a general multi-integral expression for the CF with an even number 



of spins. The case two spin operators was considered in |12] where the integration was carried 



explicitly. This result will be remembered in the Sect. |4.1| . We start with an easy example 
an proceed in order of difficulty. 



3.1 Nearest diagonal neighbors CF 

Now I review the calculation of the two-point functions. For simplicity let us consider first 
the case for nearest diagonal order-order and disorder-disorder CF Fig. 0. In term of the VO 





Cl C2 Cl C2 

Figure 3: Order-Order and Disorder-Disorder CF for nearest diagonal neighbors 



they are respectively 



E.., ^^'tr^Ns (x2^"'<l>'^'(xC2)$.(a;Ci)$'^(Ci)$<x'(C2 
g^'^gHin^ix^^''') 
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(9) 



Consider as an example (P); the reader can look at App. B in |^ for a explicit determination 
of We transform the four VO under the trace symbol in terms of fermions in the following 

way 

5^<l>.(xC2)$'^Kl)$-.(Cl)<^'"'^(C2) = 

cr 

cr 

5^^f(0)$.K2)^^^(Cl)<^>-'^Kl)<f-<x(Cl)$~'^(C2) = 

cr 

ijf{(2) / [Wil]^s<p''{wi) [Wi2] X 5^$_,K2)$-'^Kl)$-.(Cl)$~'^(C2) = 

zg^'^g'' I (f [wi2]r[wi1]ns[wi2]0^(w2)0^K), (10) 

where we changed the variables to Wi = Zi/Qf and used the notation §^ = § 2^7^, [wij] = 
Following substitutions were then made, in the second line the identity cr^ = 1, in 
the third line eqs. (|^), the intertwining relation (|^) and the second eq. in (^) in the fourth 
line and finally the unitarity conditions for the VOs. 

It only remain to take the trace of the fermion operators, which we state as a simple 
lemma, 

tr^R (x^''''(t)\w2)(t)\wi)] 

^ tr„.(.^°") (11) 

With the aid of this formula and the relation [x~^z]r[z] = — ^[zjNS we are leave with an 
integral expression for the disorder- disorder CF, 



[?i;l]Ns[wi2]Ns. (12) 

Remembering that w = exp{i7rv / 1) and ( = exp{—7ra/2I) the last expression can be written 
as: 

k f' r ■ ^ r ■ ^ 2/ /A;'\ ^ 6/(^12)) 



— av cn[v — taijcn [v — ta2) = -r 

71 J_j 7!- \k J 



TT \ k J Wia 



12) 



where H, 9, Hi and 0i (see e.g. [Q) are the Jacobi theta functions and au := ai — 02- The 
prime denotes the derivative under the argument. This integration formula appears in other 



disguise in the Appendix, see eq. ([5^ ) 
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3.2 Two-point functions with an arbitrary number of inhomo- 
geneities 

Let us call the general order-order and disorder- disorder CF by g2n{(y 1,0^2, ■■■ ,C(2n) and 
(72n(o^i! C(2, • • • 5 C(2n) respectively (when we don't want to emphasize the spectral parameters 
5^2,1 and g!^^ simply), where Q = exp(— 7raj/2/), then we have as simple generalizations of 
the trace formulas from Sect. 13.11: 



^^.f,..^, o-o-'trNs (x'^^''''F{(t, a'; fif, Ci, • • • , C2n) 

,2-Ja„ a2,..., a2n) = (g^SgK^n,,^^ (^2^NS) ^, (13) 

F(cr,o-';/ii;Ci,...,C2n) = 

<f'^'K2„)<f;.„_,K2n-l) ■ ■ ■ $'^H<2)$.Kl)$'^(Cl)<fMi(C2) ■ • • $^"-HC2„-l)$<x'(C2n), 

where the subscript in the sum above means the sum for all 1 < i < n — 1. Similarly for 
the disorder- disorder case: 

E^^trR (x2^"G(/x.;Ci,...,C2n)) 
a,,..., = (^NS,R).tr^ > (14) 

G{fii; Cl, ■ ■ ■ , C2n) = 

$Mn«2n)$'^"«2n-l) " " " K2) ^^K^^Cl) <^>-Mi (Cl) <f (C2) " " " (C2n-l) (C2„) , 

with 1 < i < n. Consider the disorder CF (the other case is similar, see |jl2[ for the 
differences). We see that if we begin in the middle of the operator expression G(/ii; Ci; • • • ; C2n) 
we have the same product of four vertex operators ([T0|) appearing before, thus we can 
substitute them by a pair of fermion operators and move them to the left by means of the 
intertwining property (^. The result is, apart from a product of elliptic functions depending 
on the spectral parameters, of the same form but with four vertex operators less as before. 
Trading in this way each time four vertex operators by two fermions operators and using the 
following antisymmetry property of the fermion algebra under the integrations, i.e., 

[wil]ji[wii]f^s [wi l][(l)^{wi),(j)^{wi)]+ = 
[wil]n[wii]m[wil]{S^{x^wi/wi) + 6^{x~^wi/wi)) = 0, 



Wl JWi 



Wl JWi 



one arrives at the expression (with the fermions ordered from left to right according to 
increasing indexes) 

(_,^NS^R)J JJ / |[^i]NsMR---[^2n-l]NsWRX J] [^^j] H </'''( ^0 • (15) 
yi<i<2n-^'^i ) l<i<j<'2n i=l 

The corresponding formula for the order-order CF is obtained by simply interchanging the 
NS with R representation and vice-versa. We continue with the analysis of this expression 
in the Section El 
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3.3 Higher even correlation functions 

By the same method of calculation of the two-point functions explained above the general- 
ization to arbitrary higher even CF is straightforward. We consider the CF with only order 
operators the purely disorder case can be obtained in the same way as before interchanging 
representations. Let us call the general 2n-CF G2n- Define first some parameters: 

• 2n spins for i between 1 and 2n at arbitrary positions in a diagonal ordered from 
left to right; 

• 2N spectral parameters Q for i between 1 and 2N, N being a positive integer; 

• 2n integers a, G 1,2, ... , 2N with Oi = 1 < 02 < 03 < ■ ■ ■ < = 2A^, such that for 
odd i the spectral parameter Q is adjacent to the right of the spin cij and when i is 
even Q is adjacent to the left of the spin 0",; 

• A set W of L := Yl^=ii^2i — Ci2i-i + 1) complex parameters Wj, note that L is an even 
integer, ordered in the following way: 

WaijWai+l, ■ ■ ■ ,'Wa2', Wa^, Wa^+l, . . . , Wa^] ... ; W^a2„_u '^a2„-i+l ) ■ ■ ■ ) ^a2„ • 

Consider the following example for a 6-point CF: In the picture we have that: n = 3, N = 7, 




Ci C2 Cs C4 Cs Ce Ct Cs Cq Cio Cii C12 C13 C14 



Figure 4: An example of a 6-point spin-spin correlation function. 



Oi = 1, a2 = 2, as = 7, 04 = 10, 05 = 13, and ag = 14. The set W of complex variables Wj 

is: Wi,W2,W7,Ws,WQ,Wio,Wi3,Wi4. 

The 2n-CF can now be written as: 



where J-'{wj} is a 'temperature' state: 



and is known, see Sect. to be a sum of products of Dirac-delta functions. 
The function Q{wj; Q} is given by the product of n blocks gj, where 

a2i-i<r<a2i a2i-i<s<a2i a2i-i<k<a2i k<j<2N 
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in the above expression the parameters r, s, k, and j are positive integers. The prime in the 
products means that we increment the parameters r and s in steps of 2, i.e. the r runs over 
odd and the s even integers respectively. 

In the example presented before we have three blocks gi, g2, and gs: 

2 14 

gl = [wil]R[w;22]Nsnn[^J^] 

i=l j<k 

10 14 

g2 = [w77]r[w88]ns[w99]r[wio10]ns JJ JJK'A;] 

j=7j<k 

g3 = [Wi313]r[Wi414]ns[W1314]. 

And g = gig2g3- 

The next step in the calculation is easy, namely to integrate the delta functions in eq. (|TB|). 
The number of integrals is reduced by a half and the final result is a clumsy sum of integrals 
expressions. Perhaps some symmetrization procedure can give more adequate results as we 
will show in the next section for the two-point functions. 



4 Explicit results for the even correlation functions 

In this section we review the general integration of the two-point functions announced in [|1^ , 
details of the calculation are relegated to the appendix. Also some examples of higher 
homogeneous nearest neighbors diagonal CF are given in the form of determinants, which 
generalize the well-known results for the diagonal two-point functions. Finally we treat the 
diagonal homogeneous four-point functions for arbitrary spins locations. 

4.1 The general two-point functions 



To integrate the expression found in Sect. p.2| we need the trace of the fermions operators 



in ([T5|) that is given by a a generalization of the lemma used before, i.e. 



tr^R(x2^^) 



with the notation 5^ (zj) = 6^{x^Wi/wj). Herepi, . . . ,p2n is some permutation of the numbers 
1,2,..., 2n, is a summation over all {2n — 1)!! permutations which satisfy the restrictions 
P2m-i < P2m foY 1 < m < Ti and P2m-i < P2m+i for 1 < m < 72 — 1 and Ep is the parity of 
the permutation. There is a similar expression for the NS case (note that the integrations 
in this case are well defined although the function 6^^{ij) has half-integral powers in the 
integration variable z) . 

Now changing the spectral parameters to their additive version i.e. interchanging Q 
with Qfj (remember that Q = exp(— 7rQ;j/2/)) and using some identities for the product of 
theta functions (see the appendix for the technicalities) we arrive at the following simple 



8 



expressions for the general two-point functions 



Pfafiian - aj)) 
92n = ^' 



Pfaffian ( v^sn 



, , , 21 fk'Y Qi'iia) , , , 2/,,,,iHi'(m) 

M«)--^(jj M«) (20) 



TT 



Expression (|T^) is symmetric in all spectral parameters and satisfies the recurrence relation 
of BaxterQ 

92n{'^l^ "2, ■ ■ ■ 5 (^2n-l, Ol2n^\ + /') = 92n-2{'^l^ ^2, " " " , a2n-2), 

+ /') = 1. 

To obtain the homogeneous correlations one has to take the limit of the expression (|l^) when 
all spectral parameters are equal which requires a generalization of the L'Hopital's rule 
for determinants already well known (see e.g. or [0), 

det|/.2i+^^-^)(0)| 
rt'^''^ = L L (2^^ 

9hom,2n det{/l(2j+2fc-3)(0)} • ^ ^ 

In the last formula /i"(0) means the n'derivative of the function h evaluated at the point 0, 
h{a) = v^sn(Q;) and /i^^o- were already defined in (pOl). The determinants have rank n and 
their elements at the row j and column k are given by h^li^a^"^^ ^\^)- the mathematical 
literature they are called persymmetric Wronskians (a terminology due to Sylvester). Be- 
cause the denominator in (|2lD can be expressed as a product one can calculate explicitly its 
homogeneous limit: 



2n-l 

n2 



det {sn(2^+'^-=')(0)} = r W (22) 



i=l 



therefore^ 



2n-l . . 

9'nL,2n = k--' n (^')-^ X det |— /.gi+^'=-3)(0)| , (23) 

This expressions are useful for studying correlations for small separations, as one only has to 
take derivatives of known elliptic functions. Several checks of eq. (^) with the results in 
were done. 



4.2 Homogeneous even-point functions for adjacent spins on a di- 
agonal 

By adjacent spins I mean nearest diagonal neighbors. I discuss the 4-point spin function and 
then the 6-point spin function the generalization afterward will be easy. 

^The integral formula conjectured by Baxter for the general two-point function can also be derived 
see H. 

^Yamada has obtained a similar expression by different methods in fl^ . 
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In the the case of the four-point CF Q is the product of two blocks gi and g2 and the 
temperature state T{wi\ is 



V W2J V «^6/ V ^5/ V ^6/ V ^^6/ V ^5 

The functions gi and g2 are 



gi = [wi1]r[u;i2][wi3][wi4][wi5][wi6] 
[w22]ns[w^23] [w24] [w25] [W26] 

g2 = [w55]r,[m;56] 
[w66]ns 

Doing the integration for the variables half of the variables Wi in ( [T^ ) gives 

J Wl J 

-i [wil]^[wi2][wi?,][wiA][wib\^s ^ i [w22]Ns[w^23][^i'24][^f25][^y26]R + 



W2 



+ i [wi1]r[wi2][wi3][wi4][wi5][wi6]r X (h [w22]ns[w23][w24][w25] 



NS- 



(24) 



Equation (|24D is the inhomogeneous 4-point function. Now I restrict to the homogeneous 
4-point function i.e., take all spectral parameters Q equal to 1. 

Let us change eq. (|2^) to the additive variables Vi {wi = exp{iTcVi/I)), but first I introduce 



the following notations: 



TC 



a^ =— dvidnvYisnvY'' for n G Z. (25) 



Note that for n < 0, a„ can be put in a different form: by the shifting v ^ v + il' in such a 
way that the poles of integrand remain outside the contour of integration and because the 
integrand is an even function the limits of integration can be leave as in (|25|) , namely 



/ dv{cnvf{snvf''''-^^ for n = l,2, 

71" J-i 



Returning to the homogeneous 4-point CF eq. (0), inserting the definitions for the different 
factors in the additive variables 



G4 = ^—Jdv{dnvY^ — ^— J dv k'^dnv cnv {snvY^ + 



where the integrals go from — / to /, the second term vanishes and with the above notation 
we have 

(^4 = flo ~ "2 a-2- (26) 
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Proceeding in the same way the 6-point CF is given by 



CLq — 2ao 0,2 0—2 ~l~ 0^2 O— 4 ~l~ 0^4 2 — ^0 ^4 0_4 



Qq (1—2 Qi—4 
02 do 0_2 



CI4 CI2 



ao 



(27) 



The general expression for G2n can be induced to be: 



G2n 



ao 
a2 



a-2 
ao 



«-2(n-l) 
fl-2(n-2) 



(28) 



0'2{n-l) 0.2{n-2) ' ' ' ^0 

Remark. The definitions of the coefficients of the above matrices is equivalent to those 
given on |Ty] pag. 199. For temperatures below the critical temperature, when < k < 1 
they are given by 



1 

2^ 



2n 



dee 



1,-1 'JF 



(29) 



deforming the contour of integration from the unit circle |^| = 1 (^ := e *^) to the real axis 
one rewrites (PD|) as 



On = - 







substituting ^ = k{snu)'^ one sees that the definition (^) coincides with (p5|). 

In a similar manner the Toeplitz determinant for the two-point function (see e.g. |ni|) 
can be derived by the method presented above, giving a strong support of the vertex operator 
method. 



4.3 Homogeneous diagonal Four-point functions 

We label the four spins from left to right as ai, (T2, as, a4 and let the distance between the 
(Ti and (72, be r and between and (T4 be s, finally let the distance between (T2 and (T3 be p, 
where the distances are measured in diagonal steps of course. Then by the same procedure 
as before, we have the following general result. 

The CF can be written as a determinant G of order r + s whose elements on the principal 
diagonal are all equal to Oq and whose elements a„ above and below the diagonal are negative 
and positive sub-indices respectively. To give its form decompose G in four matrix blocks: 



G 



B2 A2 



(30) 



Ai is the same Toeplitz matrix of order r x r for the two-point function (o"icr2) i.e. 



Ai :-- 



( ao 

a\ 



ao 



a-r+2 



\ Or-l Clr-2 " " " ^O / 
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A2 is also a Toeplitz matrix but now of order s. Bi is a rx s matrix with all its coefficients 
negative whose element in the ffist row and at the s column is equal to a-(^r+s+p)+i and with 
coefficient sub-indices growing to the left and to the bottom of the matrix 



0-{r+p)+l 



v 



-p-1 



'2-(r+p)-l ■ ■ 
■ ■ ■ 



'^-(r+s+p)+l \ 
0'-(r+s+p)+2 



Finally the matrix B2 is the transpose of but with all its coefficients positive. 

Let us give some examples for illustration. Consider the following configuration where 
each circle represents a point in the diagonal: 



• o • o o • • 

0\ O2 (T3 C"4 



The correlation function for this configuration is (with r = 2, p = 3 and s = 1): 

(35 (34 do 

An interesting example is given by the following picture 



Ox 02 



(J3 CT4 



in the figure we suppose that the distance between 02 and o"3 is n (i.e. p = n). Then the CF 
is 

O-O Q'-n-l 



these coincides with the value found in |14] using combinatorial methods 



5 Odd correlation functions 

In this section we want to give some examples of odd-point spin CF (there are not odd 
disorder CF under the critical temperature). As example of higher odd-point functions I 
treat some particular cases of the three spin CF . 

Consider ffist the three spins in adjacent diagonal positions see Fig. ^ Let us call the 3- 



Ci C2 Ca C4 
Figure 5: Adjacent 3-point correlation function. 
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point CF Ci(ai, 02, as, ^4) (the chose of the index will be clear below). Its explicit expression 
is given by 



Af 

tr„Ns fa;^^''<l>'^^(a;C4)$..«3)$'^HxC2)$..«i)$'^HCi)<^>..(C2)$'^HC3)$.3(C4) ) ,(31) 



where the normalization factor J\f is {g^g^^Yti^^x"^^^^). Now using eqs. (0) and (^) and 
starting at the middle of the operator expression in (PT|). After using the unitarity conditions 
for the VOs we arrive at 



5^a3<l>'^^(a;C4)$.3(C4) = ^7'''(id«Ns,+ - id^Ns.-) = ^^^^^o""' 



0-3 

where the unitarity relation (|^ was employed and the fermion operator i/jq^ anticommutes 
with all the fermions in the NS representation and has eigenvalue +1 in even subspace and 
— 1 at the odd subspace. Collecting all the terms we have 

Ci = -i(t i g{wi,W2) J^{wi,W2), (32) 

g(wi,W2) = [wi1]r[wi2][wi3][wi4] 
[w22]ns[w23][w24]. 



and 

tr^Ns (a;2^^V'''(w^i)0'''(^2)V'o''' 

T{wi,W2) = , -. (33) 

tr^^Ns (x^-^ j 

Evaluating the expression (|33|) gives 

T{w,,W2) = m( 5^S(x2^) + — sn {v, - V2)] , (34) 

\ W2 ITT J 

M being the magnetization^ 

Carrying out the integrations explicitly the three-point function is 

l<j<j<4 

dn (ia;24)sn (iau) [sn (m23) — cn (iai2)sn (iaia) + sn (iai2)cn (iaia)] 
■dn (ia23)sn (ia;i3)[sn (^24) — cn (iai2)sn (ia^) + sn (ia;i2)cn (20:14)] (35) 
+sn (iau) [sn (ia;34) — cn (ia;i3)sn (ia^) + sn (iQ;i3)cn (iau)] [. 



The operator represents the central spin matrix af under the Jordan- Wigner transformation, see 
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To make explicit the symmetries of this function, we employ the identities, cn {u — v) = 
cn ucn V + sn usn vdn {u — v) and dn usn. {u — v) = — cn usn v + sn ucn vdn {u — v), obtaining 



Ci = M Y\ ^ |sn (ia;i2)sn (20:34) + cn (iai3)cn (20:24) — cn (ioi4)cn (2023) + 

l<j<jr<4 

sn (20i4)sn (2023) [dn (zoi3) + dn (^024)] — sn (ioi3)sn (2024) [dn (zou) + dn (^023)] | , (36) 



eq. ( PB| ) has now symmetries under the transpositions 1 ^ 2 or 3 4 and the geometrical 
symmetry under the double transposition 1 <-> 4 and 2^3 see Fig. ^; note that the above 
expression is not invariant under the sole transposition 2 ^ 3 or 1 ^ 4. 



Taking all the Oj in (|35|)or (p^) equal and using ( ^2]) we obtain the homogeneous diagonal 
three-point function^ 



^hom (37) 

Let us investigate the case that the third spin goes to infinity the other spins remaining 
adjacent diagonal neighbors in the homogeneous case Fig. ^ By the cluster property we 
expect that the CF tends to the product of spin expectation value times the two-point 
function as the spin (T3 tends to infinity, this limit will be another check for the VOM. The 



• •00 ••■ ooo« 

(Ji (72 

Figure 6: 3-point function with one spin at arbitrary distance. 



double integral expression for the homogeneous case is then 



pn+2 rl rl 

C'^°"'.= M—^J J dvidv2dnvisn^"+hicnv2sn'^'''v2 sn {v^-V2) + -S^^ ivi-V2 +ir) ,(3J 



where n = 1 corresponds to the case of adjacent diagonal spins. The integration of the delta 
function gives 

1 2 

dvdn^v = -E{k), (39) 



where E{k) is the complete elliptic integral of second kind. The integration of the first term 
in (^) could be done in principle by residues but the most simple way seems to be to use 
the addition formula for the function sn [vi — V2), 

snt'icnt'2dnt'2 — snf2cnfidnfi 

sn [vi - V2) = 7^ — ^ ^ , 

1 — K^^sn isn ^f2 

substituting this formula we see that only the first factor is no zero on symmetry grounds 

^2n+2 sn^"'"'"^fisn^'^f2cn^f2dnfidnf2 

d^^dv2 , n_ o . _ 2.. • (40) 



TT^ J^iJ^j 1 — fc^sn 2f isn ^f2 



*This limit can also be done, of course, in the integral expression first and then integrate by residues. 
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Define the following integrals 



Jn. 



dv dnf sn ^"f, 
/ 

2n„ 



(41) 



by differentiation of the expression cn v sn v it can be proved that the integrals J„ satisfy 
the following recursion relations 

2n + 2 

In the interval — / < f i, f2 < / we can expand the denominator in (140|) 

oo 

(1 — /c^sn isn ^t>2)~"'^ = A;^-'sn^-'t>isn^-'f2, 

i=o 

putting this expression under the integrations an using the definitions of Jn one arrives to 
the following series (interchange of the series with the integrals is allowed due to the uniform 
convergence of the series in the interval of integration) 



^ ] k'^'^ Jn+j+1 {Jn+j ~ Jn+j+l] 
.j=0 



using the recursion relations and changing the index of summation one has 

2^ 27-1' 



j=n+l 



now because Jq = n hj the recursion relations 



J,- 



finally 



In= Yl 

j=n+l 



(2j-l)!! 
(2j)!! ' 

(2j - 1)! 

(2j)n 



2j-l' 



by the series definition of the elliptic integral E{k) (see e.g. 0) 



Jo = 1 - -E{k) 



TT 



we also have that 



'n ^ -'0 



(2j - 1)!! 

(2j)n 



2j-r 



therefore the final expression for the 3-point function G„ is 



(jhom = M <j 1 - ^ 



(2j - 1)!! 
(2j)!! 



2j-l 



(42) 



with n ^ oo Gn tends to Mx ^E{k) i.e. the product of the magnetization with the two-point 
function with as we expect. 
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6 Conclusion 



We have seen in the calculation of different CF how intimate is the connection of the IM 
with the theory of elliptic functions and how the vertex operator method allows to make this 
connection very directly. In this way several scattered results about the CF of the Ising model 
in the hterature could be treated in a general way. Because the method is not totally rigorous, 
specially in the definition of the CTM, several consistency checks were done: comparision 
with results already known by other methods or physical reasonable consequences. 

We have only treated purely order or disorder CF but with the same methods exposed 
in this article one can write immediately integral expressions for the mixed CF. A more 
difficult problem is to obtain a complete integration of the higher CF as was done for the 
two-point functions or alternatively to obtain possible equations satisfied by the CF like the 
ones found in H and references therein. Also the generalized Toeplitz determinants found 



in Sect. [4.2| and [4.3| could be used to see what kind of nonlinear differential equations the 



CF satisfy in the same way as was done in ||T3|. 



Acknowledgments Supported in part by DFG Sfb 288 "Differentialgeometrie und Quan- 
tenphysik" and by research funds of the Freie Universitat Berlin. 



A Integration of the two-point correlation function 

Here I give the details of the explicit integration of order-order correlation function. The 
disorder- disorder case is analogous and I will indicate in passing the modification needed. 

I begin remembering the expression obtained for the disorder-disorder CF on page ^ 
equation ([l5|) , that I write know in the following way 



fj, 

92n 



l<i<2n 



n f ks'lW X •^2^nW, (43) 



where 

P 

with the same notation as in Sect. [4.1| , and 

Q2n{Wi} ■■= (-0"[^i1]ns[w22]r ■ ■ ■ [W2n-l{2n - l)]Ns[w2n'2n]R >< H t^*-^]- 

l<«<i<2n 



As I showed in [T2|, before inserting this formula under the integral expression it is more 
convenient to antisymmetrize the factor multiplying the fermions in all the Wi variables 
(under the integrals the fermions can be anticommutated as already observed), for this we 
change the spectral parameters and integration variables to their additive version. 



'2Ik'}ii(vi-iai) , , 2Ik' QAvi - ia. 



vr Q{vi — iai) ' * V vr 0(t>j — iai) ' 
, H(t>j — iaj) 
<d>{Vi - laj) 
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The result of the antisymmetrization is 

niLi Hi(t'2i-i - ia2i-i)Qi{v2i - ia2i) ni<i<j<2n - iaj] 



ni<i<,<2n - - ^3) (Hl(O)ei(O))" 



2n 



+ perm 

2n 



Hi(o)n?:in3:i0(^. 



X 



(46) 



3) \ i=l r=l 

where by perm we mean the antisymmetrization in all vi. The same formula is valid under 
the interchange of 0i with Hi. This identity is proved showing first that both sides transform 
identically under the periods 21 and 2il' in all fj, and then through an inductive process it 
is possible show that the residues of the simple poles at e.g. vi = iai + il' are identical for 
the two sides. As a simple consequence of Liouville's theorem the rhs can only differ by an 
additive constant, which is proved to be zero by comparing the two sides at e.g. V2n = i(^2n- 
Therefore we can now write 



X 



n eK,)HK) 

l<«<i<2n 



2n 



2n 



^ ^' Hi(o) n e(t;.-.a,; 



X 



(47) 



r=l 



a,; — aj and Vij 



Vi - Vj. 



c := (ifc'Hi(0)9i(0)/7r)" is a constant factor; 

Comparing the rhs. of the above equation with the rhs. of the following classical result^ 
(see e.g. App. A in |T6[ ) 



det 

l<*<i<2n 



9i(t;,-m,)H(0)^ 
e(wi -mj)Hi(O) 



(H {[))) i<j 

Hi(0) YlQ{vi~iaj) 



we see that by multiplication of some factors dependent on the variables Vi one can write 

1 



with 



gf^{v,} = B{a,) X 
hf,{v) := 



X det 



(2n)! l<i<j<2n 

ei(i;)H(0)' 



h^{vi - aj) 



e(t;)Hi(0) 

In the order case one has instead ho-{v) := ken v. 
The function B(ai) is: 



dnv. 



(48) 
(49) 



B(a,] 



c T-r Q{ia 
(H'(0))2" ii H 



tai 



TX\fk 



sn lai 



i<j 



-1 



(50) 



Inserting eqs. (|48| ) , (|49| ) and (|50|) under the integration symbols gives (I consider now 
simultaneously the order case also) 

92n = B{ai) X S n / ^'"U 5Z n ^'^'/^(^J' ~ ^"-Pj) ^ 



^ '£:Q5(t;Ql - t;Q2 + U') ■ ■ ■ 5(t^Q(2n-l) - VQ2n + «/')• 



(51) 



^This formula is incorrectly given in reference p 
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Now each of the {2n — 1)!! terms in the sum over the deha functions gives the same result 
when integrated over their respective variables Vi, because both factors under the integrals 
are totally antisymmetric in the f^'s, therefore it is only necessary to integrate over the first 
term in the second sum i.e., 6{vi — f 2 + il')5{v^ — f 4 + W) ■ ■ ■ 5{v2n~i — V2n + H') multiplied 
by the constant (2n — l)!!/(2n)! = l/(2" n!). Integrating out the deltas one is left with (after 
renaming the integration variables) 



[Vj - mpo-+i)J 



sn lai 



i<j 



.1=1 
-1 



j=i p 



^^epha,f^{iapi - iap2) ■ ■ ■ K^^{iap{2n~i) - iap2n), 



where a = a — I' and the following definition was introduced: 



For the integration of the last formula we use (see [T 

j du cn [u + ia)cn [u + ih) = 2/[cn {ia — ib) + ik^^cn [ia — ib — il') zn {ia — ib)] (52) 

J du dn (u + za)dn (u + ib) = 2/[dn {ia — ib) + idn {ia — ib — il') zn {ia — ib)]. 



(53) 



Equations (|52|) and ( ^3]) can be shown under the condition 

-I' <a,b<I' 



for real numbers a and b. znu = Q'{u)/Q{u) is the Jacobi zeta function. Employing these 
results we have by analytic continuation, 

21 Vk, 



hrj{iai — ictj) 
h^{iai - i&j) 



vr 
21 



■ [dn ittij sn iaij — cn iaij zn iaij] + iVkcia 



[k'^cn iaij sn iaij — dn iaij zn iaij] H — ^dn 



tai 



V kiT yk 



(54) 
(55) 



The terms in square brackets in eqs. (0) and (|55|) are odd under interchanging with aj] 
the last terms are even. Inserting this expressions under the formula for only the odd 
terms survive. Now the sum in the eqs. for gfgn^ is an equivalent way of writing a pfaffian, 
thus 

-1 



92n 



HVk 



sma 



i<j 



X Pfaffian [ha,^{iai — iaj)] 



Using the equivalence of the zeta function with the logarithmic derivative of the Q{u) function 
one sees that the factors in square brackets are equivalent to the expressions on page ^ 
eqs. (^). Using now the identity 

Pfaffian (^\/ksn{iai — iaj)^ = Y\ (^y/ksn{iai — iaj)^ . 
proved by use of Liouville's theorem, see e.g. |TT[, eq. (5.8). This proves finally eq. (ITW)- 
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